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In this work we calculate the angular eigenvalues of the (n + 4)-dimensional simply rotating 
Kerr-(A)dS spheroidal harmonics using the Asymptotic Iteration Method (AIM). We make some 
'q«^" ' comparisons between this method and that of the Continued Fraction Method (CFM) and use the 

' latter to check our results. We also present analytic expressions for the small rotation limit up to 

, 0{(?) with the coefficient of each power up to 0{o?), where c = aui and oi = a^A (a is the angular 

■ velocity, uj the frequency and A the cosmological constant). 
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I. INTRODUCTION 



Recently a new method for obtaining solutions of second order ordinary differential equations (with bound poten- 
tials) has been developed called the asymptotic iteration method (AIM) [l| . The AIM provides a simple approach to 
obtainingeigenvalues of bound state problems, even for spheroidal harmonics with c a general complex number, large 
or small @, y]- It has also been shown that the AIM is closely related to the continued fractions method (CFM) [J| 
derived from an exact solution to the Schrodingcr equation via a WKB ansatz Q . A related CFM is often employed 
\l ' in numerical calculations of spheroidal eigenvalues and quasinormal modes of black hole equations Q , which is based 
^ . on the series solution method of the Hydrogen molecule ion by Jaffc (and generalised by Baber and Hasse) 0] • 

In this letter we will demonstrate that the AIM can also be applied to the generalized scalar hyper-spheroidal 
equation, Skjm{(^), derived from an (n -t- 4)-dimensional simply rotating Kerr-(A)dS angular separation equation d, [^: 



^' • -lhr,de ((1 + a cos^ 0) sin 9 cos" OdeS) + ( A.,™ - ^^^^^ - -^1^ - Kl±^ ) ^ ^ 

' sm cos" V J \ sm 6 1 + a cos^ cos"^ J 

o\'. , . . 

, where we have defined a = a A with a the angular rotation parameter. Note that A < corresponds to an asymptotic 
ILjl ■ anti-de Sitter space, whereas A > corresponds to an asymptotic de Sitter space @, and the frequency ljJ is contained 
in the dimensionless parameter c = auj. 

Higher dimensional spheroids have already been discussed by Berti et al. [l3|, who use a 3-term continued fraction 
method to solve the angular eigenvalues; however, the generalized scalar hyper-spheroida l eq uation under investigation 
there contains four regular singular points^, which leads to a 4-term recurrence relation |ll| . The simplest brute force 
approach to deal with an n-term recurrence relation is to use n Gaussian eliminations to reduce the problem to a 
tri-diagonal matrix form p^ . but this can often be very tedious in practice. 

Even in four-dimensions the Kerr-(A)dS case does not allow for a simple 3-term continued fraction relation, nev- 
ertheless an elegant method has been developed to deal with situations of this type. In this case the CFM can be 
implemented by first transforming the angular equation into the Heun form [l3|- However, such elegant techniques 
can only be applied if there are exactly four regular singular points in the equation. In contrast the AIM has an 
even broader appeal in that it can be applied somewhat^ independently of the singularity structure of the ordinary 
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^ Unlike the asymptotically flat limit (A = 0) which only has three regular singular points. 
^ Of course some manipulation is first required to put the equation into the AIM form. 
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differential equation, and thus to a larger class of equations without much modification or effort. The fact that the 
AIM may be of use when there are more than four singular points, or when other methods become prohibitively 
difficult, are compelling reasons to investigate it further. 

In this paper we use the AIM to find the eigenvalues of the simply rotating (n + 4) Kerr-(A)dS spheroidal harmonics, 
but because this equation has four regular singular points we will also take this opportunity to compare the convergence 
rate of the AIM to that of the CFM (after first transforming into the Heun form). This serves as a double check of 
our results. 

The paper is organised as follows: In section |TT] we give an overview of the AIM method and put our equation into 
the AIM form, then in section IIIII we first transform our equation into the Heun form and then describe the CFM 
method. In section |W] we analytically calculate the first three coefficients of the eigenvalue in the small c expansion, 
before finishing with some concluding remarks and analysis. 



II. THE ASYMPTOTIC ITERATION METHOD 



To write the angular equation in a form suitable for the AIM we substitute x — cos 9 and obtain: 

'n{l - x^) — x'^ 



{l-x'){l + ax^)S"{x) + 



A. 



kj7n 



+ ax{n + 2-{n + 3)x'^) - x{l + ax'^)j S'{x) 
c^{l-x'^) 771^(1 + a) jU + n-l)" 



1 + ax'^ 



1 - x^ 



Six) 



. 



(2) 



Note that the separation constant A^-j^ above corresponds to a simple eigenvalue shift in the asymptotically flat cases 
studied thus far as can be verified by setting a ~ (compare to equation (3.3) of reference [ld|)- However, if 
a such a shift is not possible, because of the non-trival 1 + ax'^ factor in the denominator. This means that the 
a — > limit will agree with the asymptotically flat case only after an eigenvalue shift. Also, note that because the 
above equation is invariant under to ^ — m we shall consider only m > 0. 

The AIM can be implemented by multiplying Skjm by the characteristic exponents (as in section IIIII for Heun's 
method); however, we have found that the most suitable form (fastest converging) is obtained by multiplying the 
angular mode function by 0, Q ■ 



Skjm{x) — (1 a; ) 2 ykjm{x) 

which leads to a differential equation in the AIM form: 

y" = Aoy' + Soy , 

where (for Kerr-(A)dS): 



Ao = - 



1 



-2 (|m| + ((|to| + 2)a;2 - l) a + l) a:^ - n (a:^ - l) [ax^ + l) 



(1 - a;2)(l + aa;2) 
1 



(3) 
(4) 

(5) 



So = 



(1 -x2)(l + aa;2) 



■■ \ + Auj,n - \m\{n + |to| + 1) - \m\ {{n + \m\ + ?,)x^ + |m|) a - 

ax^ + 1 



j{n + j - 1) 



and where the primes of y denote derivatives with respect to x. Differentiating equation ^ p times with respect to 
X, leads to: 



2/^^+^' = Xpy' + s,y , 
where the superscript p indicates the p-th derivative with respect to x and 



Ap — Ap_2 + Sp-i + AqA. 



and 



soAp_ 



For sufficiently large p the asymptotic aspect of the "method" is introduced, that is: 

Sp{x) _ Sp_l(x) _ 



\p{x) Ap_i(x) 



(3{x) , 



(6) 
(7) 

(8) 
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which leads to the general eigenfunction solution [l|: 



y{x) = exp 



(3{x')dx' 



C2 + C1 expM [Xo{x") + 2p{x")]dx"\dx' 



(9) 



for given integration constants Ci and C2, which can be determined by imposing a normalisation condition. Within 
the framework of the AIM, a sufficient condition for imposing termination of the iterations is when Sp{x) = 0, for a 
given choice of x, where Q 



dp{x) = Sp{x)Xp^i{x) - Sp^i{x)Xp{x) 



(10) 



For each value of m and k (or j), in a given (n+4)-dimensions, the roots of leads a tower of eigenvalues (m, ^i, £2, • ■ • )j 
where larger iterations give more roots and better convergence for higher i modes in the tower. 

It was noticed 0, II that the AIM converges fastest at the maximum of the potential, which in four dimensions 
occurs at a; = (even with a ^ and for general spin-s). However, in the higher dimensional case we could 
not determine the relevant Schrodinger like form and thus the maximum of the potential could not be analytically 
obtained. Nevertheless, as can be seen from the plots in Fig. [T] we found that the point a; = i = cos ^, in general, 
gave the fastest convergence. 
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FIG. 1: (Color Online) Plot of the convergence of a typical eigenvalue Ayn (71 = 1, c = 1 and a = 1) under p iterations of the 
AIM for various choices of a:: = {0.45,0.5,0.6,0.75}. Shown on the left is the eigenvalue versus p, while on the right is a log 
plot of the estimated error, \Akj,n{p) — Akjmioo)\. 
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FIG. 2: (Color Online) Plot of the convergence of a typical eigenvalue {j = m = 1, k = 7, n=l) under p iterations of the CFM. 
The eigenvalue is shown on the left, while on the right is a log plot of the estimated error, \Akjm{p) — Akjm{oo)\. 



One can also compare the relative rate of convergence between the AIM and CFM (see next section) methods for 
a typical eigenvalue by looking at the Log plots of the error, as shown in Fig[T] and Fig. [21 Due to the typically 
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exponential decrease in the error as a function of iteration number, if one is prepared to continue going to higher 
orders in the iterations of the AIM, eventually the AIM eigenvalues will exactly equal those of the CFM to any level 
of precision. 

For the purposes of consistency we calculate each eigenvalue to 10 significant figures and put in brackets the minimum 
number of iterations required to reach this precision for both methods. From Tables HI - IVl it can be seen that larger 
k modes require more iterations to achieve the same level of precision. Care should be taken when interpreting what 
information the iteration number gives us. Firstly, as one iteration of the AIM is not equivalent to one iteration of 
the CFM we can not relate this directly to computational speed. In fact although typically we need to iterate the 
AIM about twice as many times as the CFM the actual process time is often longer by factors of up to 100. That 
said we have not attempted to optimize either of the two methods here and it is not our current intention or purpose 
to do so. 

The results of the AIM for different values of to, j and k are presented in Tables [J - [V] These results shall 
be discussed in section |V] after discussing alternate approaches to obtaining the angular eigenvalue for Kerr-(A)dS 
spheroids. 

TABLE I: Selected eigenvalues, Akjm, obtained from the AIM for the Kerr-AdS case with c = 1, a = —0.05, n = 1 (extra 
dimensions) and m = 0. Numbers in brackets represent the number of iterations required to reach convergence at the quoted 
precision, where subscript A and C are shorthand for AIM and CFM respectively. 



k j=0 j=l j=2 

0.4978643318 (14)a (3)c 3.317784170 (14)a (3)c 8.201525517 (14)a (3)c 

1 8.304871188 (15)a (4)c 15.12466814 (15)a (4)c 23.93276170 (16)a (4)c 

2 23.89847347 (16)a (5)c 34.63440913 (17)a (5)c 47.35999084 (18)a (5)c 

3 47.29227791 (17)a (6)c 61.93248179 (19)a (6)c 78.56537326 (21)a (7)c 

4 78.48442957 (20)a (8)c 97.02600365 (21)a (8)c 117.5622674 (22)a (8)c 

5 117.4747381 (23)a (9)c 139.9165956 (23)a (9)c 164.3544396 (24)a (9)c 

6 164.2631560 (25)a (10)c 190.6047952 (24)a (10)c 218.9432931 (27)a (10)c 

7 218.8496664 (27)a (11)c 249.0908236 (27)a (11)c 281.3294495 (29)a (11)c 



TABLE II; Selected eigenvalues, Akjm, obtained from the AIM for the Kerr-AdS case with c = 1, a = —0.05, n = 1 (extra 

dimensions) and m = j. Numbers in brackets represent the number of iterations required to reach convergence at the quoted 
precision, where subscript A and C are shorthand for AIM and CFM respectively. 

k J = ?n = 1 j = m = 2 j' = ?n = 3 

8.302051467 (13)a (2)c 23.90067165 (13)a (2)c 47.29715054 (13)a (2)c 

1 23.89902319 (15)a (3)c 47.29444354 (15)a (3)c 78.48871605 (16)a (4)c 

2 47.29281932 (16)a (4)c 78.48633467 (18)a (5)c 117.4786244 (18)a (5)c 

3 78.48490585 (19)a (6)c 117.4764654 (15)a (6)c 164.2667830 (20)a (6)c 

4 117.4751699 (20)a (7)c 164.2647680 (22)a (7)c 218.8531204 (22)a (7)c 

5 164.2635590 (23)a (8)c 218.8512015 (24)a (8)c 281.2375960 (25)a (8)c 

6 218.8500502 (26)a (9)c 281.2357438 (26)a (9)c 351.4201872 (27)a (9)c 

7 281.2346325 (28)a (10)c 351.4183829 (29)a (10)c 429.4008807 (30)a (10)c 



III. HEUN'S METHOD FOR DE-SITTER CASE 

As we mentioned earlier we could also work with a 4-term recurrence relation directly and use Gaussian elimination 
to obtain a 3-term recurrence, which then allows for the eigenvalues to be solved using the CFM. However, if we write 
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TABLE III: Selected eigenvalues, Akjm, obtained from the AIM for the Kerr-AdS case with c — 1, a — 1, n = l (extra 
dimensions) and m — 0. Numbers in brackets represent the number of iterations, required to reach convergence at the quoted 
precision, where subscript A and C are shorthand for AIM and CFM respectively. 



k j^O j=l j=2 

0.3796279195 (32)a(7)c 3.507876364 (29)a(7)c 8.963027775 (32)a(7)c 

1 12.09034280 (29)a(8)c 21.46461733 (33)a(8)c 32.77305405 (32)a(8)c 

2 34.99526905 (35)a(9)c 50.25395559 (36)a(10)c 67.41776577 (35)a(10)c 

3 69.43716466 (39)a(11)c 90.47369703 (42)a(11)c 113.4290146 (38)a(11)c 

4 115.3704526 (40)a(12)c 142.1641442 (42)a(12)c 170.8870143 (41)a(12)c 

5 172.7901365 (44)a(13)c 205.3337886 (45)a(14)c 239.8135617 (47)a(14)c 

6 241.6949846 (48)a(15)c 279.9853860 (48)a(15)c 320.2166310 (50)a(15)c 

7 322.0845633 (51)a(15)c 366.1200572 (53)a(15)c 412.0997175 (53)a(15)c 



TABLE IV: Selected eigenvalues, A^jm, obtained from the AIM for the Kerr-AdS case with c=l, a — 1, n = l (extra 
dimensions) and m — j. Numbers in brackets represent the number of iterations, required to reach convergence at the quoted 
precision, where subscript A and C are shorthand for AIM and CFM respectively. 



k j — m = 1 j = m = 2 j = m = 3 

12.10323009 (26)a(6)c 35.37074707 (24)a(7)c 70.27222103 (25)a(7)c 

1 35.08856375 (31)a(8)c 69.80790028 (31)a(9)c 116.1867871 (28)a(8)c 

2 69.52978676 (36)a(10)c 115.7331075 (34)a(10)c 173.5944396 (32)a(10)c 

3 115.4610922 (35)a(11)c 173.1475066 (37)a(11)c 242.4916185 (37)a(12)c 

4 172.8794643 (38)a(12)c 242.0489748 (41)a(13)c 322.8761126 (38)a(13)c 

5 241.7834718 (45)a(14)c 322.4363105 (45)a(14)c 414.7467220 (45)a(15)c 

6 322.1724922 (50)a(15)c 414.3088831 (48)a(16)c 518.1027800 (50)a(16)c 

7 414.0462294 (53)a(17)c 517.6663481 (53)a(17)c 632.9438987 (53)a(18)c 



the angular equation, equation ([T]), in terms of the variable x = cos(2^^) as: 

(I - x^){2 + <5(l + x))S"{x) +(n-l~{n + 3)x + ^{l + x){n + 1 - (n + 5)x)j S'{x) 

\S{x) = 



( 



A, 



klr, 



c^{x~l) to^(1 + q;) j{j + n-l) 



2 2(2 + a(l + :r)) x-1 x + l 

and define a; = 2z — I, with the mode functions scaled by the characteristic exponents: 



Q(a-) = 2V(;2_l)V(2z)' 



viz) 



(11) 



(12) 



(note in this section we define a ^ a ^ a^A to avoid confusion with the standard Heun notation) . The angular mode 
equation can now be written in the Heun form (isl . [l5| : 
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af3z — q 



z z-1 z + lj dz z(z-l)(z + i) 



viz) = , 



where 



1 c 
7= i(2j+7i + l) , 



P=l{j + \m\+n + 3±i^) , 
^ v a 



5=1 



e=l±i- 



'a 



771 1 C C 1 

and q= — + -{j±i^){j + n + l±i^)- — 

4 4 Ja Ja 4a 



(i + |w|)(j + |m| + n + 1) - Afcj, 



(13) 

(14) 
(15) 
(16) 
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TABLE V: Selected eigenvalues, Akjm, via the AIM (we found that 32 iterations were required in every case to obtain the 
quoted precision) for different numbers of dimensions of the Kerr-dS case with c=l,a — l,m = j = k = 0. 



n 


2 


3 


4 


5 


6 


7 




0.2840487932 


0.2253670617 


0.1860599648 


0.1580675679 


0.1372036816 


0.1210960725 



with the constramt 

a + /3+l=7 + (5 + e . (17) 

Note that these results are identical to the Kerr-AdS case considered by Kodama et al. [§] by choosing a = —a?/R^ 
with c = auo. 

To compare with the AIM method we shaU use the fact that a three-term recurrence relation is guaranteed for any 
solution to Hcun's differential equation [l^, [l^ : 

aoci + /?oCo = (18) 
QfpCp+i + /3pCp + 7pCp-i =0, (p=l,2,...), (19) 



where for Kerr-(A)dS 



[p + l)(p + r - a + l)(p + r - /3 + l){p + 5) 
(2p-|-r + 2)(2p + r + 1) 
ep{p + r)(7 - ^) + [p{p + r) + a[3][2p{p + r) + i{r - 1)] 1 . ^ . 
[2p + r + l)[2p + r — 1) a 



(20) 



with 



_ {p + a-l){p + f3-l){p + -f-l){p + r-l) 

^ (2p + r-2)(2p + r-l) ' ^''> 



J + H + ^^ ■ (23) 



It may be worth mentioning that there is a removable singularity in (3p for the initial condition p = with r = 1 
(when n = 1 and j = \m\ ~ 0), which for the five-dimensional case, ti = 1, implies that this initial condition must be 
treated separately. Once a 3-term recurrence relation is obtained the eigenvalue Akjm can be found (for a given lo) 
by solving a continued fraction of the form [l3| : 

/?o-^^^..-0. (24) 

Pi— P2— P3 — 

We have used this method to compare with the AIM seen in Tables U - [V] The convergence of this method is shown 
for a typical representative eigenvalue in Fig. [21 further results are discussed in section IVl 



IV. ANALYTIC RESULTS FOR SMALL ROTATION 



It is also useful to have some analytic expressions at hand for the angular eigenvalues. For small c, these can be 
obtained by standard perturbation theory [ll| or by using eigenfunction expansion methods lid. However, a very 
convenient approach well suited to symbolic computations is the method used by Berti et al. [10|, also see reference 
[l^ . In the limit c ^ the infinite series terminates at some finite k and we are left with Q : 

Akjm - (2fc + / + |m|)(2fc + I + \m\ + n + I) . (25) 

Choosing 2k = l — {j + \m\), we find the correct c = limit: Akjm = l{l + n + l), with the constraint I > j -I- |m|. Then, 
in order to now find the small c perturbative expansion of Akjm it is convenient to use the inverted CFM, which is 
the fc*'* inversion of equation ([M)) p^ : 

r, Oik-Ilk ak-2jk~i ao7i o^k^k+i afc+i7fc+2 ,r,r-\ 
~ ~R R ■ ■ • ~W~ " It R ■ ■ ■ 

Pk-l— Pk-2~ Po Pfe+1— Pk+2 — 
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and assume a power series expansion of the form: 

p=a 

When we substitute this power series into equation (|26p the terms /„ can be found by equating powers of c (after a 
series expansion of equation ((26)) ). For the asymptotically flat case [Ifll it is very simple to go to large powers of c, 
where (in general) results to order can be obtained by going to order k = p in equation ([26]) . Unfortunately, for 
the Kerr-(A)dS case the expressions are corrmlicated by the inclusion of the curvature term a = a^A (for example the 
four dimensional case is given in reference |13|). Although results for any value of a can be stored on a computer, 
they arc too large to present on paper. Thus, given below are results up to and including 0{c^) with a further series 
expansion up to 0{a): 

/o = l{l + n + l) , (28) 

/i = — ^ ^^77- ( 2^4 + 4(n + l)P + (2m2 + 3n{n + 2) - l)f + {n + l){2m'^ +r? + 2n~- 3)1 (29) 

c [21 + n — l)(2l + n + 6) \ 



+m^(n+ l)(n + 3)-f{n^ + 2ln + An + 2l{l + 1) + 3) - .7(71- l){r? + 2ln + An + 2l{l + 1) + 3) 



c/ 



1 / (2P + 2l{n - 1) + n{n + 2) + 5) [j - I - |m|)(j - / + \m\){j + I - \m\ + n - l)(j + I + \m\ + n - 1) 



a - 



{21 + n- 3)(2Z + n- 1)^(21 + n + 1) 
c2 {2P + 2l{n + 3) + (n + 3)^) (j - ^ - |m| - 2){j - I + \m\ ~ 2)(j + / - \m\ + n + + I + \m\ + n + 1) 
^ {2l + n + 1){21 + n + 3)^{2l + n + 5) 

+ (2^ + n-l)'(2/ + n + 3) { - + - + ^ - 3) - mn + ^ + n + 1) - 21^ - Al\n + 1) 

+P (2m2 - in{n + 2) + l) - l{n + 1) {-2m^ + + 2n - 3) - 2m^{n + 3)| + 2j{j + n + I) - 2m^^ 

( 2c{-]^-jn + ] + l{l + n + l) + m^+n-l)\i 
+^ (2/ + n-l)(2^ + n + 3) ) \ ' ^'^^^ 

The singular behaviour in the denominators of /i and /2 always cancels for n = 1, because the constraint I > j + 
with 2fc = Z — (7 + |m|) being an integer, always leads to a zero in the numerator as well. This can be verified explicitly 
case by case. Note that if this expression is to be used for explicit numerical calculations then the limit must be 
chosen carefully. 

It appears that in the case of Kerr-(A)dS, odd powers of c also contribute to the spin-0 case (or tensor part of 
the graviton perturbations). Note that in the limit a — > these results do not quite agree with the results given in 
reference [l^, because of the eigenvalue shift in equation ([2]). This can simply be remedied by adding — to the /2 
term. The exact eigenvalue solution is compared with the small c expansion in Table IVTl which shows good agreement. 



TABLE VI: Comparison of the small c expansion and the exact result (via the AIM or CFM) for the Kerr-(A)dS case with 
J = 1, fc = 771 = and n = 2 (extra dimensions), for given values of c and a (results quoted to 6 s.f.). 



(c, a) 


(0.1, 0.05) 


(0.1, -0.05) 


(0.5, 0.1) 


(0.5, -0.1) 


(0.5, 0.5) 


(0.5, -0.5) 


^010 


4.01680 


3.98833 


4.09588 


4.04457 


4.17955 


3.90800 


Small c 


4.01708 


3.98863 


4.09700 


4.04585 


4.19930 


3.94355 



V. ANALYSIS & DISCUSSION 



We have calculated to ten significant figures the eigenvalues shown in Tables |T] to IIVI of the {n + 4)-dimensional 
simply rotating Kerr-(A)dS angular separation equation using the AIM. This generalises the results found in references 
0, 01 to Kerr asymptotically de-Sitter or anti de-Sitter spacetimes. Our results were also checked using the CFM. 
Although we only considered a real parameter c = aw, we could also have used a purely imaginary or complex value 
of c, and thus, the AIM may be of use for quasinormal mode analysis. For brevity we presented results for n = 1 
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extra dimensions only, but we have also checked the dependence on dimension, as can be seen in Table |V] for the 
fundamental k = mode. 

All numerics and symbolic computations used MATHEMATICA® where we found that the CFM eigenvalue so- 
lutions converged very quickly with accurate results even after a continued fraction depth of only p = 15. As a 
check of our numerics we also compared our CFM results with some independent CFM code , where we found 
identical results. (Note, that because our method is symbolic we can use the NSolve command in MATHEMATICA, 
as opposed to the method in reference [iTj that uses FindRoot). 

One point worth mentioning is that the a limit cannot be taken via Heun's method, because the recurrence 
relation (and hence the continued fraction) diverges for this case. In contrast the AIM has no such problem. The AIM 
also gives an alternative approach to obtain the eigenfunctions in terms of simple integrals, which may be useful for 
symbolic computations. Another commonly used approach would be to the use the series solution method of Leaver 
i- 

We also obtained new results for the small c expansion of the angular eigenvalue Akjm up to 0{c^^). Because of the 
complexity of these expressions, we only presented the the coefficient of each power up to O(a^), where it is interesting 
to note that the small a expansion of the small c series docs allow us to obtain the a — > limit analytically from the 
CFM. As future work it would also be interesting to compare this result with that obtained by standard perturbation 
theory [llj or eigenfunction expansion [l^ approaches. 

It is also worth mentioning that the AIM (or Heun's method) can also be applied to the case when there are two 
or more rotation parameters {non-simple) Q, where some interesting studies have already investigated spheroids in 
five-dimensions for Kerr-AdS in the near degenerate (equal rotation) limit However, each dimension must be 

considered case by case, because a general expression for (n -I- 4) dimensions has not been found. We intend to report 
on the angular spheroids (and associated radial equations) for these interesting cases in the near future. 

In conclusion, we have highlighted how the AIM can be applied to higher-dimensional scalar or tensor gravitational 
(for n > 3) spheroidal harmonics, which arise in the separation of metrics in general relativity. We have seen that 
the AIM requires very little manipulation in order to obtain a fast route to the angular eigenvalues, which may be 
useful for cases where Heun's method may not apply. However, the AIM does have some shortfalls, because although 
we did not attempt to optimise either algorithm, our implementation of the AIM was found to be much slower than 
that of the CFM. Considering that the CFM essentially involves expanding out p nested fractions, whereas the AIM 
involves taking p'^ order derivatives, this behaviour is not surprising. However, for most of the cases we considered 
only a few seconds were required to reach the desired level of accuracy and thus, timing was not a large concern. 

We hope that the AIM might be of some topical use, for example, in the angular spheroids needed in the phenomenol- 
ogy of Hawking radiation from spinning higher-dimensional black holes . We have recently used a combination of 
all the techniques discussed in this work to evaluate the angular eigenvalues, Akjrm for real c = aw, which are needed 
for the tensor graviton emission rates on a simply rotating Kerr-de Sitter black hole background in (n + 4)-dimensions 

0. 
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